Dependence of two-nucleon momentum densities on total pair momentum 
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Two-nucleon momentum distributions are calculated for the ground states of 3 He and 4 He as a 
function of the nucleons' relative and total momenta. We use variational Monte Carlo wave func- 
tions derived from a realistic Hamiltonian with two- and three-nucleon potentials. The momentum 
distribution of pp pairs is found to be much smaller than that of pn pairs for values of the relative 
momentum in the range (300 — 500) MeV/c and vanishing total momentum. However, as the total 
momentum increases to 400 MeV/c, the ratio of pp to pn pairs in this relative momentum range 
grows and approaches the limit 1/2 for 3 He and 1/4 for 4 He, corresponding to the ratio of pp to pn 
pairs in these nuclei. This behavior should be easily observable in two-nucleon knock-out processes, 
such as A(e, e'pN). 
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In a recent letter, we studied the role of tensor forces 
on the correlations between pairs of nucleons in light nu- 
clei In that work we reported calculations of the 
relative momentum distribution of pp and pn pairs with 
vanishing total momentum. We found that the strong 
spatial-spin-isospin correlations induced by the tensor 
force lead to large differences in the pp and pn distri- 
butions at moderate values of the relative momentum in 
the pair. These differences have been observed in a two- 
nucleon knockout experiment on 12 C at Jefferson Labo- 
ratory (JLab) 0. 

In this note, we report an extension of our calcula- 



tions to finite total momentum of the correlated pair for 
3 He and 4 He nuclei. This is motivated by a preliminary 
analysis of data on 3 He from the CEBAF large accep- 
tance spectrometer (CLAS) collaboration at JLab 
We find that the large differences in pp and pn distribu- 
tions gradually diminish as the center-of-mass momen- 
tum increases, until it approaches the ratio of pp to pn 
pairs for a given whole nucleus. 

The probability of finding two nucleons with relative 
momentum q and total momentum Q in isospin state 
TMt in the ground state of a nucleus is proportional to 
the density 
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where ri2 = ri — r 2 , R12 = (ri + r2)/2, and similarly 
for r' 12 and R' 12 . Here Ptm t (12) is the isospin projection 
operator, and ipjMj denotes the nuclear wave function in 
spin and spin-projection state JMj. The normalization 
is 
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where Ntm t is the number of NN pairs in state TMt- 
Obviously, integrating ptm t (<!> Q) over om y Q gives the 
probability of finding two nucleons with relative momen- 
tum q, regardless of their pair momentum Q (and vice- 
versa) . 



For this study we use variational Monte Carlo (VMC) 
wave functions, derived from a realistic Hamiltonian con- 
sisting of the Argonne vi$ two-nucleon [4| and Urbana-IX 
three-nucleon [5[ interactions (AV18/UIX). The double 
Fourier transform in Eq. ([1]) is computed by Monte Carlo 
(MC) integration. A standard Metropolis walk, guided 
by |^jMr(ri, r 2 , r 3 , . . . , r^)| 2 , is used to sample configu- 
rations [6| . For each configuration a two-dimensional grid 
of Gauss-Legendre points, Xi and Xj, is used to compute 
the Fourier transform. Instead of just moving the ip' po- 
sition (r' 12 and R' 12 ) away from a fixed ip position (ri 2 
and R12), both positions are moved symmetrically away 
from ri2 and R12, so Eq. ([!]) becomes 
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Here the polar angles of the x and X grids are also 
sampled by MC integration, with one sample per pair. 
This procedure is similar to that adopted most recently 
in studies of the 3 He(e, e'p)d and 4 He(e, e'p) 3 K reac- 
tions Q , and has the advantage of very substantially re- 
ducing the statistical errors originating from the rapidly 
oscillating nature of the integrand for large values of q 
and Q. 

The present method is computationally intensive, be- 
cause complete Gaussian integrations have to be per- 
formed for each of the configurations sampled in the ran- 
dom walk. The large range of values of x and X required 
to obtain converged results, especially for 3 He, require 
fairly large numbers of points; we used grids of up to 96 
and 80 points for x and X, respectively. We also sum 
over all pairs in a given MC sample instead of just a 
single pair. 

The pn and pp distributions at five values of the total 
momentum Q, with Q [[ q, are shown as functions of the 
relative momentum q for 3 He in Fig. [1] and for 4 He in 
Fig. [2 The statistical errors due to the MC integration 
are displayed only for the pp pairs; they are comparable 
for the pn pairs. When the total momentum vanishes, 
there is a node in the pp relative momentum distribution 
just below 2 fm , while the pn distribution has a broad 
shoulder in this region. Integration over the relative mo- 
menta in the range 1.5 — 2.5 fm -1 gives a ratio of pp to 
pn pairs R pp / pn = 0.014 ± 0.004 for 3 He, compared to 
1/2 for the whole nucleus integrated over all q and Q. 
For 4 He and Q=0 the value is R pp/pn = 0.023 ± 0.006, 
compared to 1/4 for the whole nucleus. 

The much greater magnitude of the pn momentum dis- 
tribution is due to the strong correlations induced by ten- 
sor components in the underlying NN potential. When 
Q=0, the pair and residual (A — 2) system are in a rel- 



ative S- wave. Hence, in 3 He ( 4 He), whose spin-parity is 
\ + (0 + ), pn pairs are predominantly in T—0 and 3 Si- 
3 D! (deuteron-like) states,while pp pairs are in T=l and 
1 So (quasi-bound) states |8j. The D-wave component of 
the deuteron-like pairs fills in the node in the S*-wave 
momentum distribution. 

However, for Q > 0, the two clusters may have non- 
zero orbital angular momentum and hence the pp and pn 
pairs are no longer constrained to be in the quasi-bound 
or deuteron-like states. Thus the S'-wave node in the pp 
pairs can be filled in by higher angular momentum states. 
Figures Q] and [2] show that this does happen; the node in 
the pp relative momentum distribution is rather rapidly 
filled in as Q increases. Consequently R pp / pn increases 
as Q increases, as shown in Fig. [3] for q integrated over 
1.5— 2.5 fm- 1 . 

The most direct evidence for tensor correlations in nu- 
clei comes from measurements of the deuteron structure 
functions and tensor polarization by elastic electron scat- 
tering 0. In essence, these measurements have mapped 
out the Fourier transforms of the charge densities of the 
deuteron in states with spin projections ±1 and 0, show- 
ing that they are very different. In other processes, such 
as 2 H(d, 7) 4 He 10 : ] at very low energy, or proton knock- 
out from a polarized deuteron [ll[ the effects of tensor 
correlations are more subtle and their presence is not 
easily isolated in the experimental data. This is because 
of corrections from initial or final state interactions and 
many-body terms in the transition operators. 

Some of these corrections will also affect, for instance, 
the cross sections for (e, e'pn) and (e, e'pp) knock-out pro- 
cesses in back-to-back kinematics. However, one would 
expect the contributions due to final state interactions in 
the pn and pp reactions, both between the nucleons in the 
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FIG. 1: (Color online) The pn (lines) and pp (symbols) mo- 
mentum distributions in 3 He as functions of the relative mo- 
mentum q at total pair momentum Q from to 2 fm - . 
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FIG. 2: (Color online) The pn (lines) and pp (symbols) mo- 
mentum distributions in 4 He as functions of the relative mo- 
mentum q at total pair momentum Q from to 2 fin . 
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FIG. 3: (Color online) The ratio of pp to pn pairs integrated 
over relative momentum q || Q of 1.5 to 2.5 fm -1 as a function 
of total momentum Q; red triangles are for 3 He and blue 
squares for 4 He. 

pair and between these and the nucleons in the residual 
[A — 2) system, to be of similar magnitude for relative 
momenta in the range (300 — 500) MeV/c. In particular, 
charge-exchange processes have been estimated to give 
small (< 10%) corrections to these ratios in 12 C at JLab 
kinematics |2J. Such processes, which are induced by in- 
teractions between the knocked-out pair and the residual 
cluster, could change an initial pn pair on its way out 
of the nucleus into the (detected) pp pair, thus increas- 
ing the (e,e'pp) to (e,e'pn) cross section ratios. Lastly, 
leading terms in the electromagnetic two-body current 
vanish in pp because of their isospin structure Of 
course, they will contribute in pn, but are not expected 
to produce large effects. 

The recent experiment at JLab referred to earlier has 



measured the ratio of 12 C(e, e'pn) to 12 C(e, e'pp) cross 
sections in back-to-back kinematics for relative momenta 
in the range 300—500 MeV/c to be ~ 10. These 
measurements have corroborated the results of an earlier 
analysis of a BNL experiment, which measured cross sec- 
tions for (p , pp ) and (p,ppn) processes on 12 C in similar 
kinematics j 1 31 ] . The observed enhancement in the pn to 
pp ratio is in agreement with the prediction of Ref. [l[, 
and beautifully demonstrates the crucial role that the 
tensor force plays in shaping the short-range structure of 
nuclei. 

It would be interesting to extend these measurements 
to other nuclei. In 3 He and 4 He, one would expect the 
node in the pp momentum distribution to be filled in by 
interaction effects in the final state [Tj . However, the ratio 
of pp to pn cross sections in the range (300 — 500) MeV/c 
should still reflect the dominance of the pn momentum 
distribution at these values of relative momenta. In fact, 
the analysis of JLab CLAS data on 3 He mentioned above 
suggests that this is indeed the case These data also 
seem to confirm the rapid rise of the pp to pn ratio with 
increasing total pair momentum, predicted in Fig. [3] of 
the present work. 
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